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Abstract 

Direct expressions for the magnetic anisotropy constants are given at a finite temperature from 
microscopic viewpoints. In the present derivation, it is assumed that the Hamiltonian is a linear 
function with respect to the magnetization direction. We discuss in detail the first-order con¬ 
stant K\ and show that the results reproduce previous results. We also apply our method to 
Nd 2 Fei 4 B compounds and demonstrate that the temperature dependencies of the magnetocrys¬ 
talline anisotropy constants K\, K 2 , and A 3 are successfully computed. 
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Understanding and controlling magnetic anisotropy is a central issue in developing mag¬ 
netic devices such as spin RAMs, high-density storage, and permanent magnets. The mag¬ 
netic anisotropy at room temperature (or above) is important at a practical level and should 
be evaluated in consideration of finite-temperature effects for materials whose magnetic 
properties strongly depend on temperature. This letter presents appropriate expressions for 
computing temperature-dependent magnetic anisotropy constants (MACs) and physically 
interpreting them. 

An early study of the temperature dependence of the magnetocrystalline anisotropy was 
first performed by Zener [1], which was followed by Callen and Callen [2, 3]. Considering an 
insulating ferromagnet having localized moments, they phenomenologically derived a very 
simple expression, the so-called Callen-Callen law: K n (T) = AT n (0) [M(T)/M(0)] n ^ 2n+1 ^, 
where K n (T) is the nth-order MAC at the temperature T, and M is the magnetization. This 
theory is based on the concept in which a decrease in the MAC originates from the breakdown 
of the asphericity of the electron cloud of the magnetic ions due to thermal excitation [4], 
Because a spherical electron cloud means the absence of a specific direction for the magnetic 
polarization (note that this does not mean the disappearance of the local moment), the 
temperature dependence of the MAC can be related to that of the magnetization, as shown 
above. After this pioneering work, there have been few theoretical works on the temperature 
dependence of the MACs. Much later, Skomski et al. [5-7] discussed the Callen-Callen law 
by analyzing various magnetic materials and showed that it can be applied to simple systems 
such as Fe and Co. However, they also pointed out that complex magnetic materials do not 
obey this law, such as NdCos and INcUFe^B, which exhibit a spin reorientation. Thus far, the 
Callen-Callen theory is not always satisfactory quantitatively for the range of applicability. 

In recent years, theoretical trials have been conducted to evaluate the magnetocrystalline 
anisotropy energy (MAE) at a finite temperature for spin Hamiltonians including the mag¬ 
netic anisotropy term by a Monte Carlo approach [8, 9]. Further, from a technological 
viewpoint, micro magnetic simulations of the Landau-Lifshitz-Gilbert equation including 
the random field describing thermal noise [8, 10-12] or the Landau-Lifshitz-Bloch equa¬ 
tion [13, 14] have been performed to observe the thermal stabilization of the magnetization. 

Using a first-principles approach, the MACs at T = 0 of transition-metal systems and 
rare-earth (RE) compounds have been calculated, and the results were at a satisfactory 
level compared with the experimental data within the numerical accuracy. For transition- 
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metal systems, considerable effort has been made to clarify the effects of practical factors 
such as the interfaces [15, 16] and distortion [17] in magnetic multilayers and the chemical 
disorder [18, 19] in alloys on the MAE. For the finite-temperature treatment within the 
first-principles approach, Staunton et al. [20] and Matsumoto et al. [21] have successfully 
demonstrated the temperature dependencies of the MAE of FePt and YCo 5 , respectively, 
on the basis of the functional integral method. As for RE compounds, the MACs at T = 0 
have been successfully evaluated by calculating the crystalline electric field (CEF) acting 
on the 4f electrons in the RE ions using the first-principles approach [22-24], In the 1970s, 
Wijn et al. [25] and Buschow et al. [26] calculated the magnetic stabilization energy related 
to a Sin ion in SmCos from the Helmholtz free energy using ligand field theory. Recently, 
we have theoretically studied the temperature dependence of K\ and K 2 of Nd 2 Fei 4 B and 
well-reproduced the experimental data for K t (T) using appropriate crystalline field param¬ 
eters [27]. The method employed in this work is to derive K l and K 2 numerically from the 
free energy calculated as a function of the magnetization direction. 

All approaches mentioned above have been rapidly developed in the past two decades ow¬ 
ing to the considerable progress in numerical techniques and computer performance. How¬ 
ever, it is still difficult to measure A, directly, especially at a finite temperature, because of 
the lack of an explicit quantitative expression for Ki. Given this background, we provide a 
general expression for AT* in the present work, starting from the Hamiltonian describing the 
magnetic system, and describe a method to evaluate Ki quantitatively for realistic ferro¬ 
magnetic systems. In particular, for localized electronic systems such as RE compounds or 
systems that can only be described by angular momentum operators, we demonstrate that 
the proposed method is useful and convenient for evaluating the temperature dependence of 
Ki(T). 

First, we derive the microscopic expression for the first-order MAC, AR(T), of crys¬ 
tals with ./V-fold rotational symmetry (N > 3). Introducing the Helmholtz free energy, 
F(9,(ft) := — /A 1 In ex p[— f3E n (Q, 0)] and comparing A(0, 0) with the phenomenological 
expression [28, 29], we have the relation 

K 1 (T) = ±%F(e,i,) ( 1 ) 

where /3 := 1/(/cbT) is the inverse temperature corresponding to T, and A n (0,0) is an 
energy eigenvalue of the Hamiltonian 77(0, 0). Let us consider cases where the Hamiltonian 
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is explicitly written as a function of the polar angles 9 and 0 denoting the direction of the 
magnetization relative to the c axis (the hat denotes an operator): 

U{9,<f>) := h + m(0,0) £>, (2) 

where m(6, 0) := (sin 9 cos 0, sin 9 sin 0, cos 9), and h is the angle-independent part. D is an 
operator with the transformation property of 

C' N D ± C N = e ±i %D ± , (3) 

c ] N b z c N = D\ ( 4 ) 

where D ± := D x ± i D y , and Cm is the iV-fold rotational symmetry operator around the 
symmetrical axis parallel to the c axis. The derivatives of F(9, 0) are expressed as 

d e F = Y J Z~ P{En ~ F)d eE n , (5) 

n 

d 2 e F = Y e~^~ F) [c%E n - 0 (, d e E n ) 2 ] + 0 (, %F) 2 , ( 6 ) 

n 

in terms of the derivatives of E n (9,4>). As above in the purely theoretical treatment, the 
MACs can be described in terms of E n (9,(f>) in the neighborhood of 9 — 0. Introducing 
V(9, 0) := m{9 , 0) • D — D z , we divide the Hamiltonian ? -1(9, 0) into the unperturbative part 
Hq := h + D z and the perturbative part V(6 I ,0). Here, we notice that V(0,0) —$■ 0 in the 
limit 9 — y 0 and 


C N , Ho 


= 0 . 


(7) 


To obtain A0, we perturbatively calculate E n up to the second order in V. We assume that 
the unperturbed eigenvalues e n and eigenstates | n) are given, which satisfy the unperturbed 
Schrodinger equation, Ho \n) = e n \n). Assuming that the values of e n are not degenerate, the 
perturbative energy terms are expressed as = (n|V|n), A = (n|V'P n V|n), where V n : = 
1^) (^1 /bn ~ e fc), and we have let E n (9, 0) ~ e n + A ^(9, 0) + A^ 2 ^, 0). Considering 
the symmetries in Eqs. (3) and (7), we obtain (n|.D + |n) = ( n\D + V n b + \n) = 0 such that 
d e E n (9,(j ))\ 0= o = 0 and d 2 E n (9, <f>)\ g=Q = - (n\b z \n) + [{n\b + V n b~\n) + (n\b~V n b + \n)}/2. 
Thus, from Eqs. (1), (5), and (6), we obtain 


Ki(T) 


I Y e -^ en - Fo) 

n 


(n\D z \n) + - (n\(D + V n D- + D~V n D + )\n) 


( 8 ) 
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where F 0 := F( 0, (j>). This is one of main results used to compute K\ (T) later in this letter. 
Before discussing this expression, let us further transform it. If we find an operator J 
satisfying 


[J a , ft] =ie aPl J\ (9) 

[J a 1 D^) = ie a ^D\ (10) 


we can then eliminate the matrix elements of with the help of J. Here, is the 
Levi-Civita tensor, the repeated Greek indices are summed for all cases (a = x,y,z ), the 
condition in Eq. (9) means that J is an angular momentum operator, and the condition 
in Eq. (10) geometrically means that D is a vector operator for J. Using the identity 
(n\D ± \k) = ±(e n - e k ) (n| ± (n^, h]\k), we obtain the commutator form 


Jfi(r) = 4 <[A [ A h] + h.c.]> + i Y, 

n k^n 


(n\[J + ,h\\k) | 2 + | (k\[J + , h]\n) 


€n ^-k 


(ii) 


where (• • •) denotes the statistical average in 'Hq. This is another one of the main results 
that is appropriate for describing physical pictures of the magnetic anisotropy because the 
commutators distill the essential features from the Hamiltonian. Equations (8) and (11) are 
equivalent and general for Ad- 

For example, let us apply these formulae to localized spin systems. We consider the 
single-spin Hamiltonian 


A spin := -A{S z f - 2 H ex ■ S, (12) 

where S is a spin operator; A is the single-ion MAE; H ex is the exchange held parallel to 
the magnetization direction, i.e., H cx = II ( , y rn\ and we assume the condition A < 2// ex . 
Hence, we find the correspondences D —2 H CX S and h —> —A(S Z ) 2 . Obviously, S sat¬ 
isfies the conditions in Eqs. (9) and (10) as J. Then, the commutators are calculated 
as [A+,-A(A 2 ) 2 ] = A(S Z S + + S + S z ) and [S~, [A + , -A(A 2 ) 2 ]] = 2A[S 2 -3(A 2 ) 2 ]. From 


5 



Eq. (11), 


Kf n (T) = - 


A 


M=—S 


S(S + 1) — 3 e 

( 


-P(et n ~ 


F t n ) M 2 


— Ble apl1 

* 


_ F sp in) (A 2 ( 1 - 2 Mf(S + M)(S - M + 1) 
A( 1 - 2M) - 2ET ex 


E 

M=—S 

A 2 (l + 2M) 2 (S -M)(S + M + 1) 


H(1 + 2 M) + 2H a 


( 13 ) 


where e^ in := —AM 2 — 2H ex M is an energy eigenvalue of 7Z° pin := — A(S Z ) 2 — 2H CX S Z , and 
A 0 spm := — /3 _1 In Y1m=-s ex P( — At zero temperature, we have 

Kf n {0)=As(s- 1 ^ 1 


(14) 


2) 1 + A(S — 1/2)/H ex 

It is clear that this tends to AS 2 in the classical limit S —>■ oo, as expected from Eq. (12) 
Now, we refer to the relation between i^ pm (T)/i^ pm (0) and S'(T)/S'(0). Here, the statis- 


M. The Callen- 


tically averaged spin is dehned by S(T) = ^M=_s e xp — /3(e^ m ~ -fo pm ) 

Callen law states that A'J’ pin (T)/Ji); pm (0) is [S'(T)/5'(0)] 3 . To consider this power law, we 
introduce a temperature-dependent power 

ln[i^ 1 spin (T)/A7 in (0)] 


n(T) : ln[5(T)/S(0)] 

The zero-temperature value is exactly obtained as 

6 + A(2S -3)/H c: 


a(0) = 


2 + A(2S -3)/H c _ 


(15) 


(16) 


Thus, our formula supports the Callen-Callen law: cc(0) — > 3 in the limit A/H cx —)■ 0. How¬ 
ever, we observe a deviation from the law in the classical limit because cx(0) monotonically 
decreases to 1 with respect to S. Thus, the Callen-Callen law is valid under the condi¬ 
tion AS/H ex -C 1. Figure (1) shows the temperature dependencies of the quantities with 
S = 1 and A/H cx = 0.1. We observe that ck(T) slightly deviates from the Callen-Callen law 
because of corrections from the higher-order terms with respect to A/H cx . 

The other sample is a localized 4f-electron system. Let us consider a model in which 4f 
electrons are in ligand and exchange fields [25, 28-30]: 

ACf := AZ/f • Sf + ik .cry + 2 H cx ■ Sf, (17) 
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where 7 i cry and H ex represent, respectively, the CEF and exchange field acting on the 
electrons; A is the strength of the spin-orbit interaction; and J and S are the total orbital 
and total spin angular momenta of the electrons, respectively. Assuming that H ex = H,, x m. 
the correspondences are D —$■ 2H cx St and h —> AJf •S , f+'H cry . Using J — > St, the conditions 
in Eqs. (9) and (10) are satisfied. Therefore, 


K[(T) = --(LfS 


f A 2 

t;si) + — 


E' 


-Kt'n-FS) 


E 

k^n 


(n 


Kits, 


LfSt)\k) +(n^k) 


(18) 


At first glance, we can understand that K\ vanishes in absence of the spin-orbit interaction. 
Now, we purturbatively evaluate Eq. (18) with respect to A to the second order under the 
assumption that the CEF has tetragonal symmetry. The straightforward calculation gives 


KUT) ~ — V e _/9(£ " _/) V AALZnk _ AL ™ k — ALnk , 

4 n £k ~ £n 


(19) 


where ( U cry +2H ex S f 2 ) \u n ) = e n \u n ), f := ~/3 1 In exp (-/?£„), and AL* k := | (u n \L%S?\u k ) | 2 

| (u n \Lf S^\uk) | 2 . This is a version of a localized spin system for one in itinerant electronic 
systems [31-34], Here, it is noteworthy that we can also regard Jf := L f + S f as J because 
[J“, D 1 3 ] = 0. This choice leads to another form 


K[(T) = -j ([Li, [Lt,H«,} + h.c.])' + 1 £ 

k^n 


{n\[L},U CTy \\k) | 2 + (n vy/c) 


el - el 


( 20 ) 


where we have used [Jf, XL { ■ Sf] = 0 and [Sf, 77 C ry] = 0. Thus, we can easily confirm that 
K\{T) = 0 when 77 cry = 0. The commutators can be calculated by representing 77 cry in 
terms of spherical tensor operators constructed from Jf on the basis of the Wigner-Eckart 
theorem. For simplicity, let us consider the minimum case that 


K 


cry 


= c 2 ° 


3Tf 


and 


H ex —> oo 


( 21 ) 


for A —> oo, where C 2 is a constant encoding a CEF. The commutators are calculated as 
[JjT, [J^,?7 cry ]] = —C 2 [2J 2 — 6Jf 2 ], and the second term in Eq. (20) vanishes in the limit 
of a large exchange field such that A , (T) = |C 2 (J 2 — 3Jf 2 ). Then, replacing Jf by Jf 
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on the basis of the equivalent-operator technique [35, 36] for A —* oo, we have K[(T) = 
|i ?2 {J{ — 3J f ~ 2 ). Therefore, for light RE elements, our formula reproduces 

K{(T) = -3J (J - 1/2)5°, (22) 


where B° is a CEF parameter, and J is a total angular momentum [28]. Note that in the 
limit of H ex —» 0, we observe that K[ —> 0 because the first and second terms in Eq. (20) 
cancel. 

Finally, we demonstrate numerical calculations for the temperature-dependent MACs 
(Ad, K 2 , and K 3 ) of Nd 2 Fei 4 B. Following the derivation of Eq. (8) and purturbatively calcu¬ 
lating E n {9 , 0) up to the fourth order, the microscopic expressions of K 2 and K 3 are obtained 
as 


I <2 = -,Y* e - ^ (e " -Fo) [—/f ex cy + 4 HjC® 


4 Hj(dpC® - C< 4 ») - 4ff«‘(C< 2 >C< 5 > - C< 6 >)] 


-I 


Ad 2 - J2 e ~ P{en ~ Fo) ( H ^ C L 1] - Aex 2 ^ 2 ) ) 

n 

K 3 = 2HjJ2 e ^ {en ~ Fo)C n ] i 


(23) 

(24) 


where we have defined the constants as CrP := (n\Sf\n),Cn’ := (nKSf P n Sf + S^V n S { )|n), 


( 2 ) ._ 


d 3) := (njSfVnSfln), C£> := {n\{SfVj?V n S++ S?V n S?Vjf + 2SfV n S+V n Sf + 

2 S+VJrV n S?) \n), Cl 5) := (n|(A f -iPn 2 A f + + A f + P n 2 A f -)|n), d 6) := (n| (A f + P n A+P n A f -p n A f -+ 

and ci 7) := (n\ S+V n S+V n S+V n Sf \ n). Figure 2 shows the temperature dependencies of the 
MACs for Nd 2 Fei 4 B calculated using Eqs. (8), (23), and (24). Here, we have not repeated 
the discussion of the comparison with experimental results (see Ref [27]). See Ref. [25, 28] 
for detailed computational methods for solving the eigenproblems. The expressions have 
the following advantages. (I) There are no calculation parameters. In contrast, a mesh pa¬ 
rameter is necessary for finite-difference calculations of Eq. (1). Although there are integral 
techniques using Fourier series, the spherical harmonics expansion [37], and other expansions 
in terms of various complete sets, mesh parameters are also needed. (II) The calculation 
diagonalizing the Hamiltonian matrix only needs to be performed once. If finite-difference 
or integral methods are used, the eigenvalues at the mesh points need to be computed. 


■M) 




(Ill) For materials with small MACs, high-precision results can be expected because the 
cancellation of significant digits caused by finite-difference calculations is avoided. 

In summary, we have derived direct expressions for the MACs at a finite temperature 
from microscopic viewpoints. It is only assumed that the Hamiltonian is a linear function 
with respect to the magnetization direction. We have discussed in detail the first-order 
constant, K\, and shown that the results reproduce previous results. Furthermore, we have 
successfully demonstrated the calculation of the temperature-dependent Ki, K 2l and K 3 for 
Nd 2 Fei4B. 
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T/H ex 


FIG. 1. Calculated results of <x(T) := S(T)/S( 0), K\(T) := A' 1 (T)/A' 1 (0), and the power a{T) as 
a function of the temperature T with S = 1 and A/H ex = 0.1. The Callen-Callen law predicts 
a(T) = 3. 



Temperature [K] 


FIG. 2. Calculated magnetocrystalline anisotropy constants K 2 , and K 3 as a function of 
temperature. For Nd 2 Fei 4 B, (f 2 ) = 1.001 ao 2 , (r 4 ) = 2.401 ao 4 , (r 6 ) = 12.396 ao 6 , = 295 
K/ao 2 , ^4 = —12.3 K/ao 4 , A g = —1.84 K/ao 6 , Aq = —15.9 K/ao 6 , H ex = 350 K [28], A —> 00 , and 
the lattice constants a = b = 0.881 nm, c = 1.221 nm [38]. 
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